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Abstract 



The application of the conventional saddle-point approximation to condensed 
Bose gases is thwarted by the approach of the saddle-point to the ground- 
state singularity of the grand canonical partition function. We develop and 
test a variant of the saddle-point method which takes proper care of this 
complication, and provides accurate, flexible, and computationally efficient 
access to both canonical and micro canonical statistics. Remarkably, the er- 
ror committed when naively employing the conventional approximation in 
the condensate regime turns out to be universal, that is, independent of the 
system's single-particle spectrum. The new scheme is able to cover all tem- 
peratures, including the critical temperature interval that marks the onset of 
Bose-Einstein condensation, and reveals in analytical detail how this onset 
leads to sharp features in gases with a fixed number of particles. In partic- 
ular, within the canonical ensemble the crossover from the high-temperature 
asymptotics to the condensate regime occurs in an error-function-like man- 
ner; this error function reduces to a step function when the particle number 
becomes large. Our saddle-point formulas for occupation numbers and their 
fluctuations, verified by numerical calculations, clearly bring out the special 
role played by the ground state. 
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I. INTRODUCTION 

The saddle-point method is one of the most essential tools in statistical physics ^^. 
When comparing different statistical ensembles, it is used with overwhelming success in both 
fundamental theoretical considerations and practical calculations 0,^. Yet, the conventional 
form of this usually easy-to-handle approximation fails in the case of condensed ideal Bose 
gases |§-|31; for instance, it does not yield the correct fluctuation of the number of condensate 
particles [H,^]. 

To elucidate the reason for this failure we consider an ideal Bose gas with single-particle 
energies Ei,, with z/ = 0, 1, 2, . . . labeling the individual energy eigenvalues. Since the grand 
canonical partition function H(/9,z) generates the canonical partition functions Z^if^) by 
means of the expansion 
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= Y. z'^MP) , (1) 

N=0 

each A^-particle partition function Zm{i3) can be represented, according to Cauchy's theorem, 
by a contour integral in the complex z-plane, 

where the path of integration encircles the origin counter-clockwise. As usual, (3 = l/{kBT) 
is the inverse temperature. Denoting the negative logarithm of the integrand as F{z), i.e., 
writing 
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or 
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F(z) = (A^ + 1) In ;z + ^ ln(l - ze'f^'^) , (4) 



u=0 



the saddle-point Zq is determined by the requirement that this function becomes stationary, 

dF{z 



dz 
giving 



= , (5) 

z=zo 



N + i = y —^^ . (6) 

u=Q ^0 ^ ^ 

Apart from the appearance of one extra particle on the left hand side, this is just the grand 
canonical relation between particle number A^ and fugacity zq. 



Proceeding according to folk wisdom, one then expands the logarithm F{z) quadratically 
around zq and leads the path of integration parallel to the imaginary axis over the saddle, 
relying on the fact that for large N the main contribution to the integral (0) is collected in 
the immediate vicinity of the saddle-point 0, so that the quadratic expansion should prove 
sufficient. Doing the remaining Gaussian integral, one arrives at the standard saddle-point 
approximation Z^'^ to the canonical partition functions. 



Zr'i/3) = ^ / dz exp -F-^ - ^-F^^'iz - z. 
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Here and in the following we write Z*^"-' for the n-th derivative of a function / at a saddle- 

(2) 

point; we have used F < 0. Hence, one finds 

lnZj^-P-^(/5) = --ln27r - (iV + 1) In^o - ^ ln(l - Zoe-f"'-^) 

_ilnV^ 2, (8) 

from which the canonical occupation number (11^)01 of the state a is obtained by differenti- 
ating once with respect to {—^€0), 

ainZr-) , ^nZr-^ dz, 
^ "^^'^ di-Pe^) + dz, d{-Pe^) ■ ^^^ 

As long as the gas is not condensed, the fourth term in the approximation (§) remains small 
in comparison to the third, and therefore may be neglected. Then the partial derivative 
dlnZ^^'^'' /dzo vanishes as a consequence of the saddle-point equation (|^), and one is left 
with 



(^q)c 



^r4 7 , (10) 

Z(j eP^« — 1 



so that within these approximations the canonical occupation numbers equal their grand 
canonical counterparts. 

However, this reasoning breaks down for temperatures below the onset of Bose-Einstein 
condensation, where Eq. (^ requires that z^'^^^^" — 1 be on the order of 1/A^, so that the 
third and the fourth term on the right hand side of Eq. (§) are of comparable magnitude, 
namely of the order O(lnA^). Hence, the argument that led to the familiar formula (|10|) 
becomes invalid. Moreover, inspecting the higher derivatives of F{z) at the saddle-point, 
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one finds that in the condensate regime these derivatives grow dramatically with increasing n, 
F = O^N""), casting doubt on the validity of the saddle-point approximation even if 
higher-order terms are included ^ . However, for the derivation of an asymptotic series the 
convergence properties of the formal Taylor series of F{z) are irrelevant |T^; what actually 



endangers the approximation scheme (|3) in the condensate regime is the narrow approach 
of the saddle-point to the ground-state singularity z = e^^° of the grand canonical partition 
function (0). Namely, the Gaussian integral (|^ can well represent the exact expression @ 
only if the function F{z) is free of singularities at least in those intervals where the Gaussian 
is still large, i.e., where —F (z — ZQy/2 is on the order of unity. Since F = 0{N'^), this 
observation translates into the requirement that the function F{z) should be regular at least 
in an interval of order 0{1/N) around zq. On the other hand, in the condensate regime the 
singularity at z = e^^° falls within a distance of order 0{1/N) from the saddle-point zq, again 
as a consequence of Eq. (^. It is this conflict, not the poor behavior of the Taylor expansion 
of F{z), which necessitates an approach to the contour integral (^ that is essentially more 
careful than the standard scheme (^. Interestingly, though, the dilemma does not appear 
to be overly severe — the magnitude of both conflicting intervals being of the same order 
0{1/N) — , and one might wonder already at this point whether the usual procedure can 
be saved by simple means. 

The failure of the saddle-point approximation (|^ for condensed Bose gases has been the 
subject of a long debate in the literature 0-0, with notable early contributions by Din- 



gle [11] and Fraser |12]. Various schemes have been designed for computing the number 
of condensate particles, and its fluctuation, for gases with a fixed number A^ of particles, 
without resorting to the saddle-point method. Quite recently, Navez et al. have suggested a 
statistical ensemble within which one regards the condensate as an infinite particle reservoir 
for the excited-states subsystem [Q, thus putting into shape an idea already expressed by 



Fierz |1J]. In a mathematical setting, each excited single-particle level Ea then becomes 
formally equivalent to a harmonic oscillator with frequency {Sa — £o)/h, which allows one 
to derive elegant integral representations for canonical and microcanonical expectation val- 
ues: Setting Eq = for convenience, restricting oneself to temperatures below the onset of 
condensation, and introducing the spectral Zeta function 
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with the sum running over the excited states only, the canonical number of condensate 



particles can be written as |]15 



1 rr+ioo 

K)cn = N-— dt r{t)z{f3, t)c(t) , (13) 

where C,{t) is the Riemann Zeta function. The path of integration up the complex t-plane 
lies to the right of all poles of the integrand, so that the residues of these poles, taken 
from right to left, yield a systematic expansion of the integral |jl5|. Similarly, the canonical 
mean-square fluctuation of the number of condensate particles takes the form 



1 rr+ioo 

(<^'no)c„ = — / dtT{t)Z{P, t)C{t - 1) , (14) 

and the difference between canonical and microcanonical ffuctuations is given by 

{0 no)cn - (5 ^o)mc = i cr+ioc^, .^ r, , n . ^TTZ 7^ ' V-^) 
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Yet, there are at least two reasons not to be content with this state of affairs. Firstly, since 
Eqs. ([T3|) , ([T^), and (|15D rely on the presence of a reservoir of condensate particles, they 
are blind to the onset of Bose-Einstein condensation, that is, to the sudden appearance of 
this reservoir; as a consequence of this underlying "oscillator approximation" , they do not 
allow one to discuss just how such a sharp feature can emerge in a Bose gas with a large, 
fixed number of particles. Secondly, for a given single-particle spectrum the integrals may 
not always be straightforward to evaluate; already the treatment of an anisotropic harmonic 
oscillator potential, in which case Z{(3, t) is related to Zeta functions of the Barnes type ||15|| , 



requires quite some analytical skills. Therefore, one desires a tool that works, in principle, 
for all temperatures, and is easy to use in practical calculations. 

The development of such a tool is the objective of the present paper. Following a trail 
pioneered by Dingle []T0[, we work out and test a variant of the saddle-point method that 
fulfills the above two requirements. Our guiding maxim is the same which already governed 
London's classic analysis |T^ of the condensation phenomenon: If the ground state is causing 



trouble, single it out and give it a special treatment — which, in our case, means to exempt 
the ground-state factor of the grand canonical partition function (|^) from the quadratic 
expansion performed in the usual scheme (0), and to treat that factor exactly. As will be 
demonstrated in detail, this natural strategy leads to an approach to the canonical and 
the microcanonical statistics of condensed Bose gases which is both extremely accurate and 
unsurpassed in computational ease. 

We proceed as follows: In the next section we will first explain why and how this proper 
saddle-point approximation works for the calculation of the canonical partition functions, 
and how known limiting cases are recovered, concentrating key technical details in Appen- 
dices and p. We then apply the method to computing canonical occupation numbers 
and their fluctuations, illustrating the analytical results by numerical calculations for an 



ideal Bose gas confined by an isotropic harmonic potential. In Section III we turn to the 



microcanonical ensemble, and show how the very same refined saddle-point approach allows 



one to obtain the interesting quantities almost without further effort. The final Section |V 
summarizes the most important findings. 

Although our discussion is led along the lines of the ideal Bose gas, and although some 
interesting physical insights will turn up — in particular, it will become clear why the canon- 
ical occupation numbers are well described by the expression (p^Of), and thus essentially equal 
to their grand canonical analogues, even for temperatures below the onset of condensation, 
where the previous argument had failed — , this is not primarily a work on the ideal Bose 
gas as such. Rather, major emphasis lies on the mathematics of saddle-point integrals with 
a singular integrand, of which condensed Bose gases provide perhaps the most prominent 
examples; we hope that the detailed exposition presented here will prove fruitful also in 
other areas of mathematical physics where similar problems arise. 



II. THE CANONICAL ENSEMBLE 

A. The canonical partition function 

Since the conventional approximation (^ is thwarted by the fact that the saddle-point zq, 
i.e., the solution to Eq. (H), approaches the ground-state singularity z = e^^'^ of the inte- 
grand (^ within order 0{1/N) in the condensate regime, we now exclude the ground-state 
contribution from the quadratic expansion of the logarithm d). That is, we define 

F{z) = F{z) - ln(l - ze-f^'") 

OD 

= (A^ + 1) In 2 + ^ ln(l - ze-^'^) , (16) 

u=l 

and write the canonical A^-particle partition function as 

The key idea for treating integrals of this kind, due to Dingle 0], is to let the potentially 
dangerous denominator stand as it is, and to expand only the tempered function F{z) 
around the saddle-point zq. The resulting approximation to Zi^{(3) will then be valid for all 
temperatures: For high T, when the saddle-point moves away from the singularity, it is of 
no concern whether or not the denominator is included in the Gaussian approximation; in 
the condensate regime its exclusion is crucial. 

Let us first check the behavior of the derivatives F^") . When omitting the ground-state 
term (z/ = 0) from the sum in Eq. (|ll|), the behavior of the remaining sum is governed in 
the condensate regime by the following terms (z/ = 1,2,3, .. .), each of them being about 
proportional to the n-th power of temperature. Hence, 
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.u=i v^o"'e^^— 1, 



l/n 



oc fcsT , (18) 



with approximate T-proportionality holding the better, the larger n, since large n emphasize 
the low-lying states. Now we focus on systems with single-particle energies of the form |]1^ 



£{M = ^0 + £l E ^i^t ' (19) 

where Vi = 0,1,2,... are integer quantum numbers. The dimensionless anisotropy coeffi- 
cients Cj should be of comparable magnitude, the lowest of them equaling unity. This class 
of systems contains, e.g., a gas of N ideal Bose particles confined by a c?-dimensional har- 
monic oscillator potential (s = 1), or by a (i-dimensional hard box (s = 2). For d/s > 1 and 
large N there is a sharp onset of Bose-Einstein condensation, with 

kBT oc {N,,):if (20) 

in the condensate regime; (A^cx)cn = N— (no)cn is the total number of excited particles |ir8| , p!9| . 
Thus, for largish n the Eqs. (|TB|) and (^) give 
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>gn^E(^ ^_l^,e,._ J « (A^ex)c"n/'^ for d/s>l, (21) 

resulting, by virtue of Eq. (|TID , in 

F(") = 0(A^«(")) with ^(n) =max{l,ns/4, (22) 

as fong as, besides (no)cn = 0{N), also (A^ex)cn = O(A^). To confirm this relation, the case 
n = 2 of which will be of particular interest later. Fig. |l] shows the numerically computed 
quantities 

"" = n ln(iV,.)e. ^^^^ 

(which, in the general case, should approach s/d for large n) for a gas of A^ = 10^ ideal 
Bose particles kept at temperature T = 0.5 Tg in a three-dimensional isotropic harmonic 
oscillator potential, with 



Ti-)-bLiJL\'" (24) 



denoting the condensation temperature in the large- A^-limit [0; uj is the oscillator frequency. 



As expected, r„ approaches the value 1/3 fairly rapidly with n. For comparison. Fig. 1 also 
shows the corresponding data for a gas with the same number of particles which is stored 
in a one- dimensional harmonic potential. Its temperature is T = 0.5 Tg , where 

This is a borderline case: For d = s = 1 there is no sharp onset of Bose-Einstein conden- 
sation, so that Tq merely plays the role of a characteristic temperature below which the 



ground-state population becomes significant. There are logarithmic corrections |T^ which 
keep the ratios r„ below the value s/d = 1 also for large n, but, as seen in the figure, 
even now r„ rapidly approaches a constant not too far from unity. We conclude that for 
systems of the type (p^, with d/s > 1, the 0(A^")-growth of F "" in the condensate regime 
is replaced by the somewhat milder 0(A^"''/'^)-growth of F'-"'' when going from F{z) to its 
ground-state-amputated descendant F{z), so that also the formal Taylor expansion of F{z) 
around the saddle-point zq is ill-behaved. 

But as already indicated, the properties of the Taylor series of F{z) are only of secondary 
importance. What really matters is that this function does not share the ground-state 
singularity; the singular point to be watched now is the one at z = e^^^ Since Zq < e^^°, 
the saddle-point remains separated from that singularity by at least the A^-independent 
gap e^^^ — e^^° ~ (ei — eo)/kBT. This guarantees that when quadratically expanding the 
amputated function F{z), rather than F{z), around zq, an interval of the required order 
0(l/v— F(^) = 0{N^^^'^^^'^) becomes singularity-free for sufficiently large A^; the higher the 
temperature (while remaining in the condensate regime), the smaller the gap, and the larger 
the particle number has to be. Then the Gaussian approximation to exp{—F{z)) is safe. 



As shown in detail in Appendix ^ the subsequently emerging saddle-point integral for the 
canonical partition function (|1^) can be done exactly, yielding (cf. Eq. ( |A14| ) with cr = 1) 



where 



exp(/?5o-i^(°^-l + ^^'- 
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4-P)(/3) = -^ exp{(3eo - F^'^ - 1 + ^r/^ - ^f ) D.M , (26) 



(27) 
(28) 



and D-i{rj) is a parabolic cylinder function, employing Whittaker's notation ||20|pT| . For 



discussing this unfamiliar-looking expression (EBf), which, as already remarked above, is valid 



for all temperatures, we observe [^ that D^i is related to the complementary error function 
erfc, 



D_M = exp{]f) /|erfc(-|) ; (29) 
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hence 



Zr-H/5) = exp(/3.o - F(''> - 1 + l^') ^erfc(;|) 



(30) 



For high temperatures, well above the condensation point, zq approaches zero, so that the 
parameter r], and as a consequence also f], grows without bound when A^ becomes large. 
Then we may replace the complementary error function by the leading term of its asymptotic 
expansion for large positive arguments p2| . 




^e xp(-r/V2) _ ^^^^ 



This is a special case of the approximation ( pif ) introduced in Appendix p, and implies, 
together with the further approximations ( P^ ) - (B3), that the general expression 



correctly approaches the standard saddle-point result (|^) outside the condensate phase, 

4-P-)(/?)~zJ^-P-)(/5) for high T. (32) 

This was to be expected, since when the denominator in the integrand (jT^) does not become 
small, it doesn't matter whether it is included in the quadratic approximation, as in the 
scheme (^, or treated exactly, as in the derivation of Eq. (P^). 

In the condensate regime, where e^^° — zq = 0{1/N) and F^"^^ = 0{N^^'^^) as specified 
by Eq. (^), the definition (^) of the parameter i] implies either 77 = 0{N~^^'^) or 77 = 
0{N^^'^~^), whichever is larger. Since we require d/s > 1, we conclude that in either case 
f] approaches zero for large A^, so that now r] = rj — rj^^ is a large negative number. Therefore, 
we may safely use the approximation erfc(r//-\/2) ~ 2, and arrive at 



zJ-P-)(/3)~exp(/3eo-F(°)-l 



Even when the particle number N is merely moderately large, the decrease of the parameter rj 
with temperature from large positive to large negative values can be fairly rapid, so that 



the complementary error function in the canonical A^-particle partition function (|30D acts 
as a switch, meaning that the transition from the high-temperature asymptotics (|^ to 
the condensate asymptotics (|33|) becomes quite sharp. This is confirmed by Fig. ||, which 
depicts erfc(r7/\/2) as function of temperature for a gas of A^ = 10^ ideal Bosons in a three- 
dimensional isotropic harmonic potential. In the borderline case of the one-dimensional 
harmonic potential, also indicated in the figure, there is no such sharp transition. 

It is worthwhile to discuss the condensate partition function ( ^31) a little further. Antici- 
pating that, despite the incorrect reasoning, the expression (p!o|) for the canonical occupation 
numbers will remain valid approximately even in the condensate regime, we have 

^-ie^^« = 1 + (no)c-n • (34) 

Hence, we may ehminate the saddle-point zq by setting 



f^U 



yielding 

00 -1 

z^-\m ~ eMmn.u - 1 - m.) n i . exp(/?|.o - ..d • (^«' 

The infinite product on the right hand side, describing the excited-states subsystem, equals 
the canonical partition function of a collection of infinitely many, distinguishable harmonic 
oscillators with frequencies {Si, — eo)/h (where u = 1,2,3,...), thus leading back to the 
"oscillator approximation" which has been the starting point for the derivation of the inte- 
gral representations (|I3|), (|14|), and (|I5]) in Ref. |]T5[. In contrast, the value of the present 



approximation (|36|) lies in the fact that it is not restricted to the excited states, but also 
contains the ground state explicitly. 

Taking the derivatives with respect to — /J^a, with a 7^ 0, one then finds occupation 
numbers 

(n ) -^^^^^ 

\"'a/ cr> 



1 



g/3(e„-£o) _ 1 



(37) 



and mean-square fluctuations 



2 , dHnZ^-^-^ 



(?^a)cn((^a)cn+ 1) (38) 
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of the excited states, while differentiating the logarithm of Eq. (^) with respect to — /3eo 
produces first the obvious identity 



(no)cn = N-J2 ('^-)cn (39) 

and then the important equation 

(5'no)cn = E (^'^-)cn ; (40) 



u=l 



stating that, in the condensate regime and subject to the above approximations, within the 
canonical ensemble the occupation numbers of the excited states are uncorrelated stochastic 
variables [0,0,0. 

Remarkably, the error one would have committed had one naively employed the standard 
approximation (|^ in the condensate regime, and which can be quantified only now, is not 
devastating. As explained in Appendix (cf. Eq. (|B11| )), in the low-temperature, large- 
N -regime the incorrect partition function Z^'^ (P) exceeds the correct partition function by 
merely the temperature-independent factor 1/Ri ~ 1.08444, regardless of the single-particle 
spectrum, that is, regardless of the trapping potential — implying that the error might even 
go unnoticed when carelessly taking derivatives of \nZ^'^ (P). This finding is illustrated in 
Fig. ^, again for an ideal Bose gas with A^ = 1000 particles in a three-dimensional isotropic 
harmonic potential. The figure shows the ratios of Z^'^ (P), computed according to the 
scheme (|^), of the correct approximation (p6D, and of its condensate descendant (|33D, to the 
exact canonical partition function, the latter having been obtained from the one-particle 
partition function Zi{P) by means of the familiar recursion formula 



1 ^ 
ZN{P) = j;^T.ZiikP)Z^^,{P) , Zo(/3) = l. (41) 

fc=i 

As witnessed by Fig. H, the quality of the proper approximation ( p6|) is outstanding for all 
temperatures, its low-temperature variant ( |33D performs bravely where it is expected to, 
and the standard approximation (|^) is good at high temperatures, but fails by the predicted 
factor 1/Ri in the condensate regime. For comparison. Fig. ^ shows the corresponding data 
for the one-dimensional case; here the exact A^-particle partition function is known in closed 
form |23J2^. We find features that are qualitatively similar to those in the preceding figure. 



although, with A^ = 1000, the approximation ( pq ) is not quite as good at intermediate T. 
The standard scheme (0) again is off by the same, universal factor 1/Ri at low temperatures. 

B. Canonical occupation numbers 



The derivation of the expressions (|37|) and (|39D from the condensate approximation (|33D 
to the canonical partition function serves to render that partition function plausible, but it is 
not the best one can do, since we have simply assumed the validity of Eq. (|3^). An accurate 
and fully consistent computation of these occupation numbers starts from the identity 



10 



K)en = ^^lnZ^(/5) 



1 Mdz4n- ^ ^M-pe^ 



Zn{(3) 2711 J z^ ^^ 1 — ze^]i{—j3ey) 1 — zexp(— /Se^) 
_L__L^d.exp(-G(.)). (42) 



with 



G{z) = Nhiz + Y. lii(l - ^e"^"'^) + ln(l - ze"^"") + (3eo, . (43) 

The equation that determines the saddle-point Zi for the new contour integral (H2|), namely 



°° 1 1 

^ = E ,-u..„ . + .-!.... . ^ (44) 



;2 -ig^^- - 1 zr^e^*^" 



formally looks like the grand canonical relation between particle number A^ and fugacity zi 
for a system with an extra energy level £q,. Now we have to distinguish two cases: 

If a 7^ 0, we merely have to copy the steps made in the derivation of the proper canon- 
ical partition function (^). That is, we separate the ground-state contribution from the 
exponent G{z) by defining the tempered function 



and obtain 





G{z) = 


-G{z) 


-ln(l- 


^g-/3eo 


)' 


d{- 


-%.)'^^'^- 


<f>(iz 

2m J 


exp(-G(z)) 
1 - ze-f^^o 






f^-j 


exp(/3eo 


_^(0) 


-0 



(45) 



(46) 

proceeding at once to temperatures below the onset of condensation. Then Eq. (|4^) , together 
with the previous result (Q) for Zj^{P), yields the expression 

(n«)cn = exp{F{zo) - G(^i)) {a ^ 0) (47) 

for the canonical occupation numbers of the excited states in the condensate regime. This 
accurate result is well approximated by the previous formula (pTj) : Since the level that has 
artificially been doubled differs from the ground state, we may set zq ~ -^i ~ e^^° in Eq. (^TD, 
and obtain 

(^a)cn ~ exp(ln2;o - ln(l - zqc"^^") - pe, 
1 



of3{ea-£o) — X 



(48) 

using the definitions (|I6[) and ( |45| ) of the functions F and G. 

If, however, a = 0, the exponent G{z) in Eq. (^) corresponds to a system with a doubled 
ground state. Hence, we have to temper this function accordingly, and define 
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G{z) = G{z) - ln(e^^° [1 - ze-^'"]^) . (49) 

This leads to 

d ^ ,^, If exp{-G{z) - Peo) 



-ZN{f3) = — h^ 



\2 



~-^^exp(/5.o-G(°)-2) , (50) 

where we have employed the condensate approximation ( |B9| ) to the general saddle-point 
formula ( [A 14] ), with a = 2. Thus, in the condensate regime the canonical ground-state 



occupation number is given by 

2 

(^o)cn = -=nr. 7 exp(F(2;o) - G{zi) - 1 - In 2:1) . (51) 

Zi ei^^o — 1 

For showing that this cumbersome expression actually is consistent with the familiar grand- 
canonical result in the large- A^-limit, we now have to carefully keep track of the two different 
parameters zo and zi. That is, we may set 

(^o)cn , (52) 



Zo'^ef^^o - 1 
whereas 

1 (no) 



Zl 



Ig/^^o _ 1 



(53) 



reflecting the double appearance of the ground state in Eq. (44). Hence, the argument of 
the exponential in the ground-state occupation formula ( pTD should be approximately equal 
to zero. This follows by observing 



F(^o)-G'(^i)-l-ln^i 

(iV + l)ln^-l + f:4^-^°d 

(iV + i)ln^ - 1 - fl - ^) (iV + 1 - (no)cn) , (54) 

-21 V ZqJ 



which, upon inserting 



Zl (no) en 



as obtained from Eqs. (0) and (|53D, indeed gives the required relation 

F{zo) - G{zi) -1-lnzi^O . (56) 
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For demonstrating the accuracy of the canonical formulas ( ^71) and (|5l|) , we resort once 
more to a gas of A^ = 1000 ideal Bose particles in a three-dimensional harmonic oscillator 
potential. Figure | shows the occupation number (ni)cn as a function of temperature, 
computed according to Eq. (|3), and compares these data to those that are obtained if Z]^{(3) 
and dZis[{[3)/d{—[3ei) are naively calculated from the standard saddle-point formula. The 
inset quantifies the ratios of these approximate occupation numbers to the exact ones, which 
have again been computed recursively. In the condensate regime, the naive approximation 
to both Zn{[3) and dZj^{(3)/d{—(3ei) is off by the same universal factor 1/Ri derived in 
Appendix |B|, so that this error cancels when forming their ratio: The standard saddle-point 
scheme accidentally yields the correct canonical occupation numbers of the excited states 
both above and below the onset of condensation. On the other hand, if one uses even at high 
temperatures the approximations (|33|) and (|46|), each of which is correct in the condensate 
regime only, it follows from Eq. ( pi4|) that the ratio of the individual errors is given by the 

(2) (2) 

square root of G jF . Since this ratio approaches unity in the large- A^, high-T-regime, 



Eq. ( ^71) actually is correct also at high temperatures. 

This accidental correctness is no longer met in the case of the ground-state occu- 
pation number (no)cn, depicted in Fig. 6. Now the error introduced when computing 
dZN[(3)/d{—(3eQ) at low temperatures within the standard scheme is given by the factor 
l/i?2 (see Appendix ^, so that the resulting value of (no)cn is too small by the con- 
stant factor R1/R2 ~ 0.96106 in the condensate regime, whereas Eq. ( pTD yields the cor- 
rect data. Outside the condensate regime the standard approximation becomes correct, 
whereas Eq. ( |B14| ) reveals that the condensate approximation (|5TD is wrong by the factor 
M2/M1 ~ i?2/(-Ri"\/2) ~ 0.73576. Needless to say, if one deduces the canonical occupation 
numbers directly from the ratios of the saddle-point approximations ( |A14| ) to the respective 



contour integrals, instead of invoking their high- and low-temperature limits (pq ) and (|B9[) , 
one obtains expressions that are valid for all temperatures. 

C. Canonical fluctuations 

The calculation of the canonical mean-square fluctuations {S'^na)cn of the occupation 
numbers now directly parallels that of the occupation numbers themselves: Starting from 
the identity 

2 \ d^ In Zn 



2 



5(-/3£a) 
1 dZ^ ( 1 dZ^y 2 1 r ( -rri ^\ f... 

we are left with the task to evaluate the further contour integral 



where 



H(z) = (A^ - 1) In z + ^ ln(l - ze''^'") + 2 ln(l - ze"^^") + 2/?£„ . (59) 



u=0 
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-i./d.exp(-^(.))=-i./dz.^^P(-^(^)) 



As a consequence of the second derivative performed in Eq. (^71), the state with energy Ea 
now has formally been tripled, giving the equation 

oo 1 2 

N -l = y —^^ + —^^ (60) 

for the new saddle-point Z2- Thus, for a 7^ the tempered version of the function H{z) 
becomes 

H{z) = H{z) - ln(l - ze-^'°) , (61) 

implying 



27ri I V ^ V 2TTi I 1 — ze~f^^" 

-exp(/3£o-i^^°^-l) (62) 

in the condensate regime; hence 

{S^n^)cn = (^a)cn " ML + 2 exp{F{zo) - H{z2)) (a ^ 0) . (63) 

Setting zq ~ -^2 ~ 6^*^°, similar to the reasoning underlying Eq. (^8]) , it can be seen that the 
third term on the right hand side is close to 2{na)l^, so that this expression (^) properly 
reduces to the famihar Eq. (|38|). 

If a = 0, we disentangle the three-fold ground-state contribution from H{z) by defining 

H{z) = H{z) - ln(e2'^""[l - ze-^^'^f) , (64) 

and get 

^ Uz exp(-^(.)) = ^ Uz'-^^^^tlM^lpl 

~K^*;;^)''=''K*« -*'"'- 3). (65) 

employing the condensate approximation ( |B9| ) with cr = 3. Therefore, the canonical fluctu- 
ation of the number of condensate particles now takes the form 

(5'no)cn = (no)cn " (no)e'„ + —^^ ^ exp(F(^o) " H{z2) " 2) . (66) 

We check the results (|^) and (^) again for A^ = 1000 ideal Bosons in an isotropic 
harmonic oscillator potential. Figure |^ shows the root-mean-square fluctuation {6ni)cn = 
{5'^ni)l[^ as obtained from Eq. (p3|), and from the standard saddle-point approximations to 



the three individual terms on the right hand side of the identity (|57|) . For each term we 
have the same accidental correctness of the standard saddle-point result in the condensate 
regime, and of the condensate approximation at high temperatures, as already described 
for the occupation numbers of the excited states; the inset, which depicts the ratios of 

14 



the two approximations to the recursively calculated exact fluctuation, confirms that either 
approximation is correct at both high and low temperatures. 

In the case of the condensate fiuctuation, however, the situation is quite different. As 
witnessed by the inset in Fig. H, our formula (|66|) gives the correct fluctuation of the number 
of ground-state particles in the condensate regime. Since the canonical mean-square fluctu- 
ation ((5^no)cn has to vanish for zero temperature, when all A^ particles occupy the ground 
state, the third term on the right hand side of Eq. ( |66D approaches A^^ — A^ for T — > 0. 
When naively using the standard saddle-point scheme, the results for the individual terms 
in Eq. (|57D are incorrect by factors R1/R2, {Ri/R2)'^, and R1/R3, respectively, with the 
universal numbers R^ determined in Appendix |B[ Hence, the naive scheme yields spurious 
mean-square ground-state fluctuations 



{S'no) 



spur 
en 



[RjRs - {R,/R2f]N' + [R1/R2 - Ri/R3]N 
0.02438 A^^ + 0.01304 A^ 



(67) 



for T — >■ 0, giving, for instance, the incorrect r.m.s.-fluctuation (5%)cn"'^ ~ 156 for A^ = 1000, 
in precise agreement with what is observed in Fig. |^. In the opposite regime, that is, for 
high temperatures, the standard scheme becomes correct. Then, since (no) en ^ 1, the 
condensate approximation (^) is off by roughly the same factor R2/{Ri\/2) ^ 0.73576 that 
also determines the error of (no)cn itself. 

It is now also illuminating to compare the saddle-point scheme developed here to the 
oscillator approximation that has led in Ref. [15] to the integral representations ([T3|) - 
(P!3|). This latter approximation cannot cope with the Bose-Einstein transition, that is, 
its validity remains restricted to the condensate regime, since it rests on the flction of an 
inflnite reservoir of condensate particles [|l3l. This very feature, however, is what allows one 
to derive closed expressions for the condensate fluctuations, provided the pole structure of 
the spectral Zeta function ([T2|) in the complex t-plane is known. For an ideal Bose gas in 
a three-dimensional isotropic harmonic potential, this function can be written in terms of 
Riemann Zeta functions, 



Z{p,t) = {(3huo 



i-t 



\^^' 



2) + 2C(t 



i) + C(t) 



(68) 



having shifted the ground-state energy to Eq = 0. Taking into account the three rightmost 
poles of the product r(t)Z{f3,t)({t — 1), located at t = 3, t = 2, and t = 1, Eq. ( P^ then 
gives 



(5'%)c 



+ 



+ 




r 



4+C(2) 



(69) 



for ksT /{huj) ^ 1, where 7 ~ 0.57722 is Euler's constant. In contrast, the saddle-point 
approach requires the numerical determination of the three saddle-point parameters zq, 
Zi, and Z2 from Eqs. (pf), (H), and (pO)? respectively, but this effort is rewarded by the 
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possibility to monitor the fluctuations for all temperatures, including the transition regime. 
As detailed in Fig. |^, sufficiently below the onset of condensation both the saddle-point 
approximation (^6]) and the oscillator approximation (^) yield excellent agreement with 
the exact condensate fluctuation, even for particle numbers as low as A^ = 1000. 

III. THE MICROCANONICAL ENSEMBLE 

A. The microcanonical partition function 

For extending the techniques developed in the previous section to the microcanonical 
ensemble, we write the grand canonical partition function (|l]) as 

OO -| oo 

where the microcanonical partition functions Q{E, N) denote the number of microstates 
accessible to a thermally isolated A^-particle system with total excitation energy E. Setting 
the ground-state energy equal to zero, we now assume that all single-particle energies e^, 
and hence also all possible excitation energies, can be represented reasonably well as integer 
multiples of a basic quantum hu. Introducing the variable 

X = e"^^" (71) 

and writing E /{Tioj) = m, Eq. (^) takes the form of a double power series, 

oo -| oo oo 

^Eix,z), (72) 

with g^y indicating the degree of degeneracy of the single-particle energy level uhuj, and 
Q{m,N) = Q{mhuj,N). With the help of two suitable contours which encircle the origin 
of the complex x- and z-plane, respectively, the desired partition functions can be isolated 
from this series by means of the identity 

Iff -j^ CO -|^ 

= -n^—yl y da; y dz exp (-F(x, z)) , (73) 

where 

oo 

F(x,z) = (m + l)lna; + (A^ + l)lnz + ^c/^ln(l-;zx^) . (74) 



The saddle-point (xq, zq) for the double contour integral ([73| ) now is determined by the 
simultaneous solution of the two equations 
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dF{x, z) 



dx 



x=xo,z=zo 



dF{x, z) 
dz 



(75) 



X=XQ,Z=Z(i 



reading explicitly 



m + 1 = ^ Qr^—T-^ 

i/=0 ^0 ^0 -■- 



^ + i = E^^-T-i 



v=0 ^0 ^0 



1 



(76) 



As long as the gas is not condensed, we may safely use the standard approximation 
scheme [§. That is, we may expand the function (|7^ quadratically around the saddle- 
point. 



F{x, z) ^ F(xo, zo) + \f^'''\x - XoY + F^'^'\x - xo){z - z^) + \f^'''\z - z,f 



(77) 



with 



pi.r,s) ^ d^+'F{x,z) 



dx^dz^ 

substitute x — xq = iw and z — zq = iu, and get 

1 



x=xo,z=zo 



rt{m,N) 



{2ny 



exp 



(^-F{xo,Zo] 



X r^dw r^du e^vC-F^'-'^w' + F^'''^wu+-F^'''\' 



1 



exp 



{-n 



27iVD 
where D is the functional determinant 



Xq, Zq 



(78) 



D = det I 



-p(2,0) -^(1,1) 
-^(1,1) -p(0,2) 



-^(2,0)-^{0,2) _ /-p{l,l) 



(79) 



For temperatures below the onset of condensation, however, this procedure is invalid, because 
then the second of the saddle-point equations (|76D dictates that Zq differs from the ground- 
state singularity z = 1 of the grand canonical partition function ([72| ) merely by a quantity 
of order 0{1/N), exactly as in the canonical case. Hence, we have to proceed according to 
the insight accumulated there for computing Q{m, N) in the condensate regime, and have 
to exempt the ground-state factor from the the quadratic expansion: Defining 



F{x,z) = F{x,z) -ln(l ~ z) 
we have to start from the representation 



n{m,N) 



(27ri) 



1 I A I A exp{-F{x,z)) 

w^nn — r^z — 



(80) 



(81) 
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The expansion of the tempered function F{x, z) around (xq, zq) then yields 



F(x, z) ^ F{xo, zo) + -F(2'0)(x - xo) 



+ 



1 



where we have used 



F(o,i) + ^(1.1) (^ _ ^^)j (^ _ ^^) + iF(°'2)(^ - Zof 



F(i'O) = 



(82) 



(83) 



It is now of key interest to compare the magnitude of the two terms in the square brackets 



of Eq. 



which constitute the coefficient of (z 



Zo) 



On the one hand, we have 



oz 



= 20 



1 - Zo ' 



(84) 



hence F^^'^' = 0{N) in the condensate regime. On the other hand, relevant contributions to 
the integral over x are collected from an interval of order Of 1/v— F(^^) around xq. Thus, 
the relevant F^^'^\x — a;o) are on the order of F^^'^y \^^F^^ , with 



F(1'1) = 


-Y.9u 




v=l 




oo 


^(2,0) ^ 


-T.9u 



z/x^ ^ 



ZqXqj 



v2 ' 



V ZqXq 



v-1 



v=\ 



ZQxix 



(85) 



If we focus again on systems of the type ( [T9| ) — i.e., if both the quantum fiuj and the 
weights g^, have been adjusted accordingly — , we may repeat the reasoning that has led 
to the canonical estimate (^Tj), and conclude that both F^^'^) and F^^'^^ are of the order 
(9(^j\r2s/d^_ Therefore, F^^'^)(x — xq) = 0{N^^'^) for relevant x. To be honest, this can be 
taken as a rather crude guideline only, in the same sense as the quantities r„ displayed 
in Fig. |l] have not yet approached the expected value 1/3 for n = 2. In fact, numerical 
calculations for the three-dimensional isotropic harmonic oscillator potential reveal that for 
reasonably large A^ the expression F^^'^^V—F^^^^^ is about proportional to (A^ex)mc^, the 
square root of the total number of excited particles in a microcanonical setting, instead of 
being proportional to (A^ex)mc^- Nonetheless, the above estimate indicates that for d/s > 1 
and large A^ we may neglect F^^'^\x — xo) against F^^'^\ This implies a drastic simplification 
of the analysis, because then the remaining saddle-point integral factorizes: Leading both 
contours parallel to the respective imaginary axis over the saddle, we are left with 



n{m,N) 



exp{-F{xo,zo)) f^o+ 



27ii 

1 rZQ+ioo 

X-— / dz 

ZTli J ZQ—ioo 



/ dxexp --F(''°)(x-xo)' 

exp(-F(0'i)(z - Zo) - |F(°'')(^ - zof 



\-z 



(86) 



The first of these integrals is standard, the second is precisely of the type worked out in 
Appendices ^ and ^. Thus, without any further labor we obtain 



nK^)^ --'-^^^g)-^ (87) 

the saddle-point approximation to the microcanonical partition functions in the condensate 
regime. 

B. Microcanonical occupation numbers and their fluctuations 

The computation of microcanonical occupation numbers {na)mc in the condensate regime 
is a matter of routine by now, so we merely need to sketch the main steps. Denoting the 
number of microstates of an isolated iV-particle system with total excitation energy mhcj 
and with exactly n^ particles occupying a given single-particle state with energy ahu as 
Ta{na; rn, N), we have 

N 

J2 ^a{n„; m, N) = n{m, N) (88) 

and 

^^"^-^ f^KAO ■ ^^^^ 

Introducing the symbol d/d{x°'), where the overbar is meant to indicate that the partial 
derivative acts on only one of the ga-^oXd degenerate states with energy ahu, the first 
microcanonical moments are generated from the grand canonical partition function ( [72[ ) by 
means of the identity 

J) oo oo I N \ 

o;"^^ s(x, ^) = E ^"^ E ^"^ I E ^«r,K; ^, ^) 1 • (so) 



Af=0 m=0 \na=0 



Hence, 



N Iff 

E naTa{na;m,N) = jdx jd 



^(x,z) 



71^=0 



(27ri)2 / / x'^+iz^+i d{x^) 



= 7^— T^ydx^dz exp(-G(x, z) , (91) 



where 



G(a;,z) = (m + l-a)lnx + A^lnz + E^i^ln(l-zx'')+ln(l-zx") . (92) 

The evaluation of the integral ( |9TD first requires the knowledge of its saddle-point (zi, zi), 
obtained by simultaneously solving the two equations 

^ y a 

m 



+ 1 - « = E 9u -1 -y r + ,1 _» r 

i/=0 ^1 ^1 ~ J- Zi Xi — i 



°° 1 1 
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If then a 7^ 0, we define the tempered function 

G{x,z) = G{x,z)-\n{l-z) , (94) 



yielding 



^ r t AT\ ^ Ia I a exp{-G{x,z)) 

2^ ra„r„(ra„;m,A/) = fdx fdz 

nt^o (27rz)2/ / 1-z 

exp{-G{xi,zi) - 1) 



V-27rG(2-0) 
in direct analogy to Eq. (|87D, with 



(95) 



z/^^riXi ^ a^ZiXi ^ 



^1 (1 - 2;ix^)^ (1 - zix^y 
If, however, a = 0, we have to account for ground-state doubling. In this case we define 

G{x,z)=G{x,z)-2\n{l-z) (97) 

and invoke Eq. ( P9|) with cr = 2, resulting in 

v^ ^ , ^^, 1 f,f, exp(-G(x,z)) 

j:^noTo{no;m,N) = j—^^ fdx fdz ^^ _ ^^, 

(98) 



2 exp(-G(xi,Zi) -2) 



1-^1 V-27rG(2.o) 



To give at least one application of these formulas. Fig. IC shows the microcanonical ground- 
state occupation number (no)mc as a function of the microcanonical temperature for A^ = 
1000 ideal Bose particles in an isotropic harmonic oscillator potential, as computed from 
Eqs. ( pTl) and (|98|) according to Eq. (^91). The microcanonical temperature does not differ 
significantly from the canonical one |P; the inset quantifies the ratio of canonical to grand 
canonical, and of microcanonical to grand canonical occupation numbers. As expected, the 
relative differences between the ground-state occupation numbers in the three ensembles are 
on the order of 1/A^. 

For calculating the corresponding microcanonical mean-square fluctuations {S'^na)mc 
within the saddle-point approximation, we exploit the identity 

oo oo I N \ 

= E ^"^ E ^" E nlT^{n^- m, N) (99) 

Ar=0 m=0 \nQ=0 / 

which immediately leads to the analogue of Eq. ( pTj) , namely 

('^^na)mc = (na)mc " (^a)mc + ^/ ^x ^^t^-^^ j^^jdz exp(^-H{x, z)) . (100) 
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The newly appearing integral is defined by 

■^-^ <fdx <fdz exp[-H{x,z) 

giving 

OD 

H{x,z) = {m + l-2a)lnx + {N -l)\nz + J2gJn{l-zx'') + 2\n{l-zx'') . (102) 

u=0 

Hence, its saddle-point (x2, -22) is found by simultaneously solving the two equations 

m+l-2a = }^ Qu—r-^ — - + — r— ^ 7 ' 

i/=o ^2 2:2 — i ^2 2:2 — i 



1 

~Q Z2 'X2 — i Z2 ^X2 ' 



N-l = J2g, „, _, , + -1 -. , • (103) 

,, n '^9 -^9 -L ^^9 Xo 1 



The usual distinction follows: If a 7^ 0, extracting the ground-state contribution from the 
function ( |102|) means introducing 



resulting in 



H{x,z) = H{x,z)-\n{l~z) , (104) 



1 i^ /• exp(-/j'(x,z)) exp(-ij'(x2, 2:2) - 1) 



with 



^^-223^2 ^ o «^-22a;2 ^ 



^«N2 ■ 



(106) 



rr(2,0) _ _ y^ '^ ^2J-2 o "^ ^2-^-2 

.^/'(1-^2X^)' (l-;^2a;f; 
If a = 0, we define instead 

i7(x,2) =;ff(x,;z)-31n(l-z) , (107) 

and invoke Eq. ( [B9| ) once more, now for cr = 3, to arrive at 

^ /d. /d. ^^P^^^^'^^^ ^ 1 ^ exp(-i7(x2,.2)-3) _ g^ 

(27rz)2/ 7 (1-^)3 2(1-Z2)2 V-27r/7(2-0) ^ ^ 



Collecting the results (|83), (H), and (^) for a ^ 0, or (0), (H]), and (|T08|) for a = 0, 
Eq. ( |100| ) then allows one to determine the fluctuations. An example for such a calculation 
is depicted in Fig. O: The heavy solid line is the root-mean-square fluctuation ((5r2o)mc = 



(5^?T,o)mc ^s obtained from the above saddle-point scheme for a gas of 10^ ideal Bosons in the 
usual isotropic oscillator potential. For comparison, when evaluating for the same system 
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the formula 
one finds 



((5^^o)cn - ((^^^0)1 



T5|) for the difference ((5^no)cn — ((^^^0)1110 up to terms of the order ksT /{Tiu) 
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2 C(4) 2C(3) 



(109) 



This, together with Eq. (|5UD, yields a closed expression for the microcanonical condensate 
fluctuations of an ideal Bose gas in a three-dimensional isotropic harmonic trap in the 
oscillator approximation, that is, under the assumption of an infinite reservoir of condensate 
particles: 



{6^no)r 



0.64366 



kuT 



huj 



+ 
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ksT 



+ 1.85443 



knT' 



hto 



0.55531 In 



ksT' 



0.96969 



ksT 



(110) 



The dashed line in Fig. [Tl| corresponds to the leading term of this approximation [|I3| , whereas 
the thin line also takes the next-to-leading order into account. Even for 10^ particles, finite- 
size effects are still visible in the condensate fluctuation; the leading-order term of Eq. ( |110| ) 
alone yields only modest agreement with the saddle-point result. After accounting for the 
dominant corrections, the agreement becomes close to perfect: Even on the scale of the 
inset, the saddle-point result is indistinguishable from the oscillator approximation ( p.lO| ). 



IV. DISCUSSION 



The necessity to abandon the usual saddle-point scheme when exploring canonical or 
microcanonical statistics of condensed Bose gases with A^ particles is brought about by a 
characteristic dilemma. On the one hand, the approach of the saddle-point to the ground- 
state singularity at z = e^'^° of the grand canonical partition function within order 0(1/A^) 
may be taken as the very hallmark of Bose-Einstein condensation; on the other, the custom- 
ary Gaussian approximation requires that intervals of order 0{1/N) around the saddle-point 
stay clear of singularities. The solution to this problem almost suggests itself: If one exempts 
the ground-state factor of the grand canonical partition function from the Gaussian expan- 
sion and treats that factor exactly, but proceeds as usual otherwise, then the singular point 
that now decides the fate of the approximation is the one produced by the first excited state 
a.t z = e^^i. Since the saddle-point remains pinned below e^^°, it remains separated from 



the decisive singularity a.t z = e^^^ by an A^-independent gap. This gap is wide enough to 
get the approximation going if the particle number is sufficiently large, because the required 
interval of regularity shrinks with increasing A^. The representation (|17|) of the canonical 
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partition functions can be viewed as the prototype integral expressing this strategy; the 
other canonical and microcanonical quantities computed in this work constitute nothing but 
variations of the same mathematical theme. 

The success of this amended saddle-point method hinges on the fact that the emerging 
integrals with singular integrands can be done exactly; as explained in Appendix ^ they 
lead directly to parabolic cylinder functions. Thus, we have accomplished the two goals set in 
the Introduction: The results ( |A14| ) are easy to use, and provide bona fide approximations 
to partition functions, occupation numbers, and their fluctuations which are valid at all 
temperatures — not only in the high-temperature limit or in the condensate regime, but 
also in the critical temperature range that witnesses the onset of condensation. In particular, 
the sharpness of the drop of the complementary error function contained in the canonical 
partition function (|30D allows one to precisely assess the sharpness of this onset in a Bose 
gas with a finite, fixed number of particles. 

The observation that the interval of regularity claimed by the Gaussian expansion around 
the saddle-point is of the same order 0{1/N) as the distance of the saddle-point from the 
ground-state singularity, meaning that the original conflict is not too large, reflects itself in 
the noteworthy fact that the error of the conventional scheme in the condensate regime is 
merely a temperature- and system-independent multiplicative constant. The discovery, made 
in Appendix P, that this constant approaches unity when the multiplicity of the ground- 
state pole is increased fits in nicely: Increasing that multiplicity amounts to considering a 
Bose gas with a multiple-degenerate ground state and thus drives the saddle-point away from 
the ground-state singularity — the occupation number of each individual of the degenerate 
states is lowered — , thereby lessening the error of the naive approach. 

It is also of interest to compare the workload implied by the proper saddle-point method 
to that required by other techniques aiming at canonical or microcanonical statistics. Exact 
recursion relations like Eq. (|4l|) or their microcanonical analogues are invaluable for treating 
relatively small samples with not substantially more than about 1000 particles |T^,5^, but 



the computation of, e.g., microcanonical condensate fluctuations by such means for a gas 
with 10^ Bosons, as presented in Fig. |rT|, is entirely out of the question. The integral 
representations (p!3D - (|1^), on the other hand, immediately yield analytical expressions 
for condensate occupation numbers and fluctuations, provided the pole structure of the 
spectral Zeta function (0) is known, but they do not allow one to monitor the system at 
the onset of condensation. In contrast, saddle-point techniques always require one numerical 
step — finding the saddle-point as the root of the respective saddle-point equation — , but 
once this has been done, the formalism yields all quantities of interest, by means of the 
ever-same formulas, without further hardship. Therefore, we may conclude that despite all 
reservations [§,0 piled up in more than half a century since Schubert's incisive comments 0, 
it really is the saddle-point method which, if executed properly, provides the most powerful 
approach to the statistical mechanics of isolated, condensed ideal Bose gases. 

Having an instrument that reliable and flexible at one's disposal is certainly not merely 
of mathematical value, but may also have some bearing on experiments with Bose-Einstein 
condensates of dilute atomic vapors which are now becoming routine. These experiments 
are mainly done in isolated harmonic traps, in the regime of vapor densities where the 
atomic interactions, quantified by the s-wave scattering length a, can be considered as weak: 
Denoting the atom mass as m and the characteristic trap frequency as cu, and defining the 
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oscillator length L = Jli/{muj) which quantifies the extension of the trap's ground state, 
one has N{a/L)^ <^ 1 under typical conditions. For example, a = 5.4 nm for *^Rb p9 



giving N{a/ Vf = 8 • 10~^ for a sample of iV = 10^ condensate atoms in a shallow trap with 
uj = 100 s~^, while Na/L = 2 ■ 10^. These two relations place the system in the Bogoliubov 
regime, traditionally associated with the notion of a weakly interacting Bose gas. However, 
it is feasible to prepare even more weakly interacting samples, either by tuning the scattering 
length with the help of external magnetic fields |^ , or by working with spin-polarized atomic 
hydrogen |^], which features the unusually low triplet scattering length a = 0.0648 nm 



Thus, a/L = 2.6 ■ 10~® in a harmonic trap with uj = 100 s~^, so that even for A^ ^ 400 000 
hydrogen atoms one finds Na/L ^ 1, indicating the crossover regime from the ideal to 
the Bogoliubov gas. This crossover should manifest itself, in a non-trivial manner, in the 



behavior of the condensate fiuctuations |^ , which also determine what one may aptly term 



the minimum linewidth of an atom laser" [B1,B^. It would therefore be of substantial 



importance to study condensate fiuctuations of very weakly interacting Bose gases, that is, 
of systems intermediate between the ideal gas and the Bogoliubov gas, and to probe whether 
the difference ( p.09|) between the canonical and the microcanonical ensemble remains visible 



there; in general, this difference should show a pronounced dependence on the trap type |]T5 
In this way, an old, apparently purely academic issue — the non-equivalence of statistical 
ensembles in the condensate regime — suddenly pops up at the forefront of topical research, 
in the theory of the atom laser. Seen from the experimental angle, such an enterprise is on 
the verge of becoming possible; on the theoretical side, the first requirement is a tool for 
routinely computing ideal condensate fiuctuations within the different ensembles, for traps 
with various geometries. This tool is available now. 

APPENDIX A: ACCURATE SADDLE-POINT APPROXIMATIONS FOR 

BOSE- TYPE INTEGRALS 

In Section |n| we have met contour integrals of the form 

with positive integer a, and a saddle-point lying too close to the singularity at z = e^^° 
for the standard approximation (|^) to be viable. In this appendix we derive the proper 
saddle-point approximation to these integrals, following a suggestion by Dingle |l^ . 



Writing, in accordance with our previous notation, the negative logarithm of the full 
integrand as f{z), 

Jiz) = f{z) + ia- l)f3so + a ln(l - ^e"^^") , (A2) 

the saddle-point z^, is determined by the equation 



dz 



. (A3) 



In the large- A^-limit, this equation corresponds to the grand canonical relation between 
particle number N and fugacity z^, for a system with a a-fold degenerate ground state. 
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Substituting z = 2* + m, so that the saddle-point is found at m = 0, and writing the 
difference between the singular point and the saddle-point as 



we have 



and 



Mo = e^'" - z, , (A4) 



/(z* + u) = f{z^ + u) - /3eo + a\n{uo - u) (A5) 



e*°Jj,exp(-/(,-+.))^ (A6) 



2Tii J [uq — uY 

Defining /(") = f^''\z^), Eqs. (^ and (^ immediately yield 



(1) = 


-a -— \n{uo - 
du 


-u) 




a 






Uq 





u=0 



(A7) 



moreover, we require Z*-^-' <^ 0. Expanding /(z* + u) up to second order around the saddle- 
point — with a first derivative f^^^ which does not, as in the conventional approximation (^, 
vanish, but instead becomes large when uq is small, as in the condensate regime — , then 
leading the path of integration over the saddle at m = 0, we obtain the approximation 



2vrU(-P-) = (-i)-exp(/?£o-/ 



(0) 



X / ''^du expf-- - i/(2)n^ ) {u - u,r , (A8) 



which, upon substituting 

V 

u = i , + Mo , (A9) 

becomes 

0--1 



27r.eP-) = (-l)-exp(/3.o - /^^^ - a + ^v') (^^ 

/+oo+iri / 1 o\ 
dv expiiriv v 1 v^'' , (AlO) 



with 



r/^MoV^^ (All) 



a 



r] = ri . (A12) 

V 
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The integral occurring here is closely related to Whittaker's parabolic cylinder function 
D^^ijj), namely |36| 



f dv exp(ir]v - -vA v'" = V2^r''e~i'^^ D_^(ji) , (A13) 

where C runs from — cxd to +00, passing v = from above. Thus, the saddle-point approxi- 
mation to the integrals ( [AID finally takes the form 



cr-l 



X exp(/?eo - /^°) -a+^-v'- \f) D_,{fj) . (A14) 



APPENDIX B: HIGH AND LOW-TEMPERATURE APPROXIMATIONS FOR 

BOSE-TYPE INTEGRALS 



The approximation (|A14|) is valid regardless whether or not the saddle-point z^, lies close 



to e^^°, that is, for all temperatures. Therefore, it should adopt a more simple form in the 
condensate regime, where z~^e^^° — 1 is of the order 0{a/N), and should merge into the 
standard saddle-point formula in the high-temperature limit, where z* approaches zero. 

In the case of high temperatures, the saddle-point moves away from the ground-state 
singularity; the distance Uq defined in Eq. (|A4|) approaches unity, Uq ^ 1. Hence, both 
parameters i] and rj = r] — cr/rj introduced in Eqs. ( |A11|) and (|A12| ) adopt large positive 



values. Then the asymptotic expansion of the parabolic cylinder functions gives ||20 



D.M ~ '-^^^^^ ; (Bi) 



^<T 



moreover, we have 



1 .2 1 2 



2^' - -r/^ ~ a . (B2) 



Observing that, as a consequence of uq — > 1, 



77^r7~^-/(2), (B3) 

/(") ^ /(o) _ peo , (B4) 

7^'^ = f^'^ - ^uf ~ /(^) , (B5) 



we find 






exp (-7^°^) 
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(B6) 



This is just the expected result: For high temperatures, i.e., when z^ stays sufficiently far 
away from e^^", we recover the formula provided by the usual approximation scheme (|^). 

In the opposite limit, that is, in the condensate regime, we infer y — /*^^^ = 0(iV^*-^^''^) 
from Eq. (^), whereas Uq is of order 0{a/N). Since C(2)/2 = max{l/2, s/rf} for systems 
of the type (pJ]), and we have required s/d < 1, we find that Uq goes to zero faster than 
.J—f{2) increases when A^ becomes large. Hence, 77 = uqJ —f^'^'> approaches zero for large A^. 
This, in turn, implies that the argument rj of the parabolic cylinder functions now is a large 
negative number, so that pO | 



Since 



D^M - 7^^i-Vr~' exp(f /4) . (B7) 

(ct-1)! 



fj ^==, (B8) 



the low-temperature limit of the approximation ([A 141) becomes 



^^■'-7;7^Q°"-p(*»-/'°'-). (B9) 



which has been used heavily in Sees. |I| and |Tl|. 

It is of interest to recast this expression (B^) into a form which lends itself to a direct 



comparison with the now incorrect standard saddle-point formula. To this end we exploit 
that, since Mq — > 0, the second derivative of the full function / at the saddle-point will be 
dominated by the singular ground-state contribution, giving 

r ~ -4 ■ (BIO) 

With this additional approximation, and utilizing the identity / = /^^^ — (3eo + alnuo, 
Eq. (^ yields 

/^■^-T^-oexp(-/°^-a 
((T-1)! V 

^"-'e-- exp(-/^ 

~ — - v Zna — - 

- i?. ^^^tz!i) . (Bii) 

This is a most intriguing observation: In the large-A^, low-temperature limit, i.e., in the 
condensate regime, the result of the properly executed saddle-point approximation to the 
Bose-type integral ( |A1| ) differs from the standard saddle-point formula by a temperature- 
and system-independent factor i?^-, namely 



R, = V2^-, ttt; (B12) 

(a-lj! 



some numerical values of R^ are listed in the following table. 
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a 


Ra 


1 


0.92214 


2 


0.95950 


3 


0.97270 


5 


0.98349 


10 


0.99170 



Recalling Stirling's formula for (cr — 1)!, one immediately realizes that these universal 
renormalization factors R^ approach unity when the singularity index a is increased. This 
finding, which might appear paradoxical at first, has a simple explanation: In a system with 
a cr-fold degenerate ground state, each individual of these states takes only l/a-th of the 
particles that a non-degenerate ground state would have to carry. Therefore, z~^e^^^^ — 1 is of 
the order 0{a/N): The larger a, the farther away is the saddle-point z^, from the singularity 
at z = e^"^", and the better is the standard procedure. 

On the other hand, in the high-temperature regime, where uq ~^ 1, the left hand side of 
Eq. ([B9D can be written as 



d'^-ie"'^ 



(a-1)! 



exp -/ 



f(o) 



M, 



exp(-/°^^ 



-2./^) 



(B13) 



Thus, when using the low-temperature approximation (p9| ) in the high-temperature regime, 
the result is incorrect by the factor 



M^ = R^ 



rio . 



(B14) 
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This integral representation can be obtained from Eq. (19.5.4) in Ref. [Q, observing 
that the parabolic cylinder function U{a,z) equals the Whittaker function D_a-i/2{z). 
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FIGURES 



FIG. 1. Circles: ratios r„, as defined by Eqs. (|2l| ) and (|23|), for a gas of A^ = 10 ideal Bose par- 

(3) 
tides with temperature T = 0.5Tq which is confined by an isotropic three-dimensional harmonic 

oscillator potential. According to Eq. (|2l|) , one expects r„ ~ 1/3 for largish n. The diamonds 

indicate the corresponding data for a gas with same number of particles kept at T = 0.5 Tq in a 

one-dimensional harmonic trap. 



FIG. 2. Complementary error function eiic{r] / ^/2) , with the temperature-dependent parameter 
rj defined by Eqs. (27) and (28), for a gas of A^ = 10^ ideal Bose particles in a i/iree-dimensional 



(3) 

isotropic harmonic oscillator potential (full line; here the reference temperature Tq equals Tq ' as 
given by Eq. (p4|)), and in a one-dimensional harmonic potential (dashed line; with Tq = Tg as 
in Eq. (^)). When erfc(r7/v2) approaches zero, the standard high-temperature result (^ holds; 
when it approaches two, Eq. (|33| ) provides the correct partition function. 



FIG. 3. Ratio of the standard saddle-point result (^) (solid line approaching the value 
1/Ri ~ 1.08444 at low temperatures), of the proper approximation (p^) (solid line everywhere 
close to unity), and of its low-temperature descendant ( |33| ) (dashed line), to the exact canonical 
partition function Zn{(3), for a gas of A^ = 1000 ideal Bose particles in a three-dimensional iso- 
tropic harmonic potential. The reference temperature Tq is given by Eq. (^). Note the impressive 
overall performance of the approximation (p^). 

FIG. 4. As Fig. 0, now for a gas of A^ = 1000 ideal Bose particles in a one-dimensional harmonic 
potential. The reference temperature Tq is given by Eq. (p5|). Note that for T/Tq -^ the ratio of 
the standard approximation (^ and the exact data approaches the same value 1/Ri ~ 1.08444 as 
met in the three-dimensional case. 



FIG. 5. Canonical occupation number (ni)cn as obtained from Eq. ( |47| ) (full line), and from the 
ratio ofdZN{P)/d{—f3ei) to Ziy{(3), both computed from the standard saddle-point formula (dashed 
line), for A'^ = 1000 ideal Bose particles in an isotropic three-dimensional harmonic potential. The 
inset shows the respective ratios of these approximate occupation numbers to the exact ones. The 
energy level ei is three-fold degenerate; the data shown here correspond to an individual state. In 
this and all following figures, the reference temperature Tq is given by Eq. (p^). 



FIG. 6. Canonical occupation number (no)cn as obtained from Eq. ( |5lD (full line), and from 
the ratio of dZiy{(3)/d{—(3eQ) to Ziy{f3), both computed from the standard saddle-point formula 
(dashed line), for the same system as considered in Fig. |5|. The inset shows the respective ratios of 
these approximate occupation numbers to the exact ones. The error of the standard approximation 
in the condensate regime is determined by the ratio R1/R2 ~ 0.96106; the error of the condensate 
approximation ( pi] ) in the high-temperature regime is given by the factor i?2/(-Riv2) ~ 0.73576. 
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FIG. 7. Canonical r.m.s. -fluctuation {6ni)cn as obtained from Eq. (33) (full line), and from the 
standard saddle-point approximations to the individual terms in Eq. ( |57| ) (dashed line), for the 
same system as considered in Fig. |5[ The inset shows the respective ratios of these approximate 
fluctuations to the recursively computed exact one. As in Fig. |5|, the data refer to only one of the 
three states with energy ei. 



FIG. 8. Canonical r.m.s. -fluctuation (5no)cn as obtained from Eq. ( pq) (full line), and from the 
standard saddle-point approximations to the individual terms in Eq. ( |57| ) (dashed line), for the 
same system as considered in Fig. ^ The inset shows the respective ratios of these approximate 
fluctuations to the recursively computed exact one. The error of the standard scheme for T — > is 
determined by Eq. (|67|), whereas the error of the condensate approximation for high temperatures 
is given by the factor (i?2/(^i ^2))^/^ « 0.85776. 

FIG. 9. Canonical r.m.s-fluctuation {5no)cn for the same system as considered in Fig. g. The 
exact, recursively computed fluctuation (long-dashed line) is compared to the data obtained from 
the saddle-point formula (|6^) (heavy full line), to the approximation provided by the leading, 
r^-proportional term in Eq. ( p9| ) (short-dashed line), and to the prediction made by the oscillator 
approximation (^) with terms up to order T (thin line). The inset emphasizes the outstanding 
accuracy of both the proper saddle-point method and the oscillator approximation. 

FIG. 10. Micro canonical ground-state occupation number (no)mc for a gas of A'^ = 1000 ideal 
Bosons in a three-dimensional isotropic harmonic oscillator potential, as computed from the sad- 
dle-point approximations ( p7| ) and (|98[). The inset shows the ratio (no)cn/(f^o)gc of the canonical 
to the grand canonical occupation numbers (short dashes; with exact, recursively computed canon- 
ical data), and the ratio (?^o)mc/(^o)gc of the microcanonical to the grand canonical values (long 
dashes) . 

FIG. 11. Microcanonical r.m.s. -condensate fluctuation (5no)mc for a gas of A^ = 10^ ideal 
Bosons in a three-dimensional isotropic harmonic potential, as computed from the saddle-point 
approximation (full line). The dashed line corresponds to only the leading term of the oscillator 
approximation, cf. Eq. ( |110D ; the thin line, visible only in the upper right corner, also takes the 
next-to-leading term into account. Even in the inset, the result of the saddle-point calculation is 
indistinguishable from this oscillator approximation. 
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